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Abstract 

As a rough model for the collective motions of cells and organisms we develop here the 
statistical mechanics of swarms of self-propelled particles. Our approach is closely related 
to the recently developed theory of active Brownian motion and the theory of canonical- 
dissipative systems. Free motion and motion of a swarms confined in an external field is 
studied. Briefly the case of particles confined on a ring and interacting by repulsive forces is 
studied. In more detail we investigate self-confinement by Morse-type attracting forces. We 
begin with pairs N = 2; the attractors and distribution functions are discussed, then the case 
A' > 2 is discussed. Simulations for several dynamical modes of swarms of active Brownian 
particles interacting by Morse forces are presented. In particular we study rotations, drift, 
fluctuations of shape and cluster formation. 

1 Introduction 

We consider here the collective modes and the distribution functions of finite systems of free 
particles, of particles confined in an external field and systems of interacting Brownian particles 
including active friction. This is considered as a rough model for the collective motion of swarms 
of cells and organisms PElEl^. 

First we will discuss the basic model of active Brownian particles and introduce several models 
of active friction. The energy supply is modeled by a velocity-dependent function. Beside the 
classical model of negative friction due to Rayleigh [SBHlEj we will discuss the model of Schienbein 
and Gruler fS] and the so-called depot model which has been derived from concrete assumptions 
about the energy supply pi llO[[TT|ll2lll3j . We include into the model also a weak global coupling 
to the velocity of the center of mass. 

Then the motion of a swarm of noninteracting particles which is confined in an external field 
is studied. Furthermore the case of particles confined on a ring and interacting by repulsive forces 
is investigated. Next we will investigate self-confinement by Morse-type attracting forces. The 
Morse potential is defined by 



- 1 , (1) 



where r is the distance of two interacting particles. This simple model of interactions describes 
repulsive and attractive interactions similar to the Lennard-Jones potential. The Morse potential 
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which has the right physical shape is however much more useful with respect to analytical treat- 
ment. This is well-known from quantum mechanics, where even exact solutions of the Schrodinger 
equation for the Morse potential are known. Since we are very much interested in analytical so- 
lutions we decided to use the Morse potential also as an appropriate model for the description of 
the essentially non-physical interactions in swarms. 

We begin with the case of pairs N — 2 with Morse interactions, the attractors of motion and 
the distribution functions are discussed in detail. New solutions for the dynamical modes are 
presented in the deterministic and in the stochastic description. The border of stability of the 
modes is discussed. Finally the investigation is extended to systems with a large number of active 
Brownian particles interacting by Morse forces. In particular we will study the left/right rotations 
of pairs, clusters and swarms. We will show that the collective motion of large clusters of driven 
Brownian particles reminds very much the typical modes of collective motions in swarms of living 
entities. 

In the theoretical part we will make use of methods developed in the context of the theory of 
the so-called canonical dissipative systems [nill5lll6Tll7lll8j . This theory is in close relation to our 
recently developed theory of active Brownian particles P?l ll0lll9llllll2()j . 



2 Dissipative forces and equations of motion 

Let us consider two-dimensonal systems of N point masses m with the numbers 1, 2, i, ...N. We 
assume that the masses m are connected by pair interactions. If the distance between the mass i 
and the mass j is denoted by rij = ri — rj, the force is —U'{r). 

The dynamics of the system is given by the following equations of motion 

Here the dissipative forces are expressed in the form 

F,{v,) = -mj{v^^)v, (3) 

The function 7 denotes a velocity-dependent friction, which possibly has a negative part. The 
second term on the r.h.s. models a weak global coupling to the average velocity of the swarm 

i 

This way the dynamics of our Brownian particles is determined by Langevin equations with dissi- 
pative contributions. The Langevin equations contain as usually a stochastic force with strength 
D and a (5-correlated time dependence. 

{m)=0; (C.W^,(0)=^(^-i')% (5) 

In the case of thermal equilibrium systems we have ^{v) — jo — const.. In the general case where 
the friction is velocity dependent we will assume that the friction is monotonically increasing with 
the velocity and converges to 70 at large velocities. In the following we will use the following 
ansatz based on the depot model for the energy supply p iTUll^ 

where c, d, q are certain positive constants characterizing the energy flows from the depot to the 
particle These assumptions lead to the dissipative force law 
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Dependent on the parameters 70, c, d, and q the dissipative force function may have one zero at 
V = or two more zeros with 

vl = ^C, where — -1 (8) 

is a bifurcation parameter. In the case C > a finite characteristic velocity vq exists. Then we 
speak about active particles. We see that for \v\ < vq, i.e. in the range of small velocities the 
dissipative force is positive, i.e. the particle is provided with additional free energy. Hence, slow 
particles are accelerated, while the motion of fast particles is damped. In the case of small values 
of the parameter ( we get the well-known law of Rayleigh (see jITj ') 

F^mjo(^C~^ V (9) 

and in the opposite case of large values of C we get the law derived by Schienbein and Gruler from 
experiments for the dynamics of cells 

F = m7o (^^^^l^v (10) 

This way, we have shown that our depot model covers several interesting limiting cases discussed 
earlier in the literature 



3 Brownian particles in parabolic confinement 

Let us consider a finite many-particle system consisting of N particles which are driven by dissi- 
pative forces and are confined by external forces. Interactions and global coupling are neglected 
so far. It is important to understand first the motion of non-interacting Brownian particles before 
going to interacting systems. For the Langevin equation of motion we postulate 

mv, + WU = F,{vi) + V2D^^{t). (11) 

Here U is the potential of the external forces and Fi is the velocity-dependent dissipative force 
discussed above. Sometimes we will use later units with m = 1. We begin with the case of a 
parabolic external potential 

U = U{r) = ^colr' (12) 

For non-interacting particles the A'^— particle distribution factorizes and it is sufficient to treat 
the one-particle problem. The one-particle distribution function is at the same time the relative 
density of a swarm of non-interacting particles moving in the external field. The Fokker-Planck 
equation and the standard methods to obtain the stationary distribution functions are in detail 
discussed by Klimontovich jH| for the one-dimensional case. These methods were extended to the 
case of two-dimensional oscillators in |llj . In the mentioned work however only approximative 
solutions could be presented. Here we make use of the fact that our system is of canonical- 
dissipative type and will present a new solution which is exact at small noise Let us first 

consider the deterministic motion. Under stationary conditions the particles have to obey the 
requirement of balance between centrifugal and attracting forces jr = |[/'(r)|. For the harmonic 
potential this leads to the stationary radius r-g = vq/loq. Actually the particles are moving in the 
neighborhood of two limit cycle orbits which have the projections given above and are located on 
two surfaces in the four-dimensional space corresponding to the angular momenta L — ±Vq/ujo- 
This way the probability is concentrated on two closed curves in the four-dimensional phase space 
which are similar to tires In order to find the solution we go step by step and use several 
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relations valid on the limit cycle orbit. For example we have \L\ = H/loq. Further the potential 
and the kinetic energy are exactly equal 

Using this relation we may replace 7(1"^) by j(H/m), this is exact on the I.e. and a good approxi- 
mation near to it. This way we obtain a solution for the probability density which is concentrated 



flixi,X2,Vi,V2) = C 





2D 


■ H ' 




exp 




2c 


2kT 



(14) 



This expression however is not yet the correct solution since the corresponding angular momentum 
is not perpendicular to the coordinate plane. In order to concentrate the probability of the two 
tires which are observed in the simulations we make use of the relation H = ±LuJo also valid 
on the I.e.. Here 

L ^ m {xiV2 - X2V1) (15) 

is the angular momentum. On the limit cycles it has the values L — zbmrowo- In the stochastic 
case the angular momentum is distributed around L ~ ±H /ljq. We assume that this distribution 
is Gaussian. Combining now both parts of the distribution in such a way that the dissipative 
probability flow disappears, we obtain the following approximation for the stationary distribution 
function 

fo{xi,X2,vi,V2) = fi{ex.p[~ai{H - uJoL)]+exp[-ai{H + ujqL)]} (16) 

In the case — s- this approximative solution is concentrated around the two limit cycles. From 
the condition that the mean energy remains unshifted we get 

"1 = ^-^0 (17 

Summarizing this results we may conclude that the shape of the probability in the four-dimensional 
phase space is well understood from theory and numerical simulations. However the available 
approximate solutions reflect only the limit of small noise so far. Exact solutions for finite noise 
level are not yet known. 

Before we proceed to the general case of interacting particles let us make a study of a limiting 
case of interacting systems which may be treated in an elementary way. In the case of very strong 
driving, i.e. Vq 3> kT all Brownian particles will move on large orbits tq — vq/luq. This way the 
particles form a ring. For more general confining forces with radial symmetry U (r) the radius ro 
is determined by the condition of equilibrium between centrifugal and centripetal forces 

^ = \U'iro)\ (18) 

The effect of weak noise will be an equal distribution on the ring and some small dispersion 
perpendicular to it. Switching on arbitrary repulsive forces between the Brownian particles, e.g. 
an exponential repulsion (Toda forces) between neighboring particles leads to a kind of lattice on 
the ring with the average distance 

'-^ («) 

This way we come to the interesting conclusion that in the limit of strong driving, confined particles 
with repulsive interactions will always form rotating ring lattices. 

A more detailed study of the case of Toda rings has been given in our earlier work |13lll2l05) . 
Replacing the repelling exponential Toda forces by Morse forces which have an attracting tail 
beside rotating ring lattices also clustering phenomena on rings are observed 
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4 Dynamics of self-confined clusters of driven particles 



The study of the many-body dynamics of interacting driven particles is an extremely difficult task. 
Therefore we begin our investigation with the study of 2 driven Brownian particles |26| which are 
self-confined by attracting forces. In this relatively simple case we observe already several basic 
features of the dynamics of swarms. Let us consider two Brownian particles which are pairwise 
bound by a radial pair potential U{ri — ^2), as a concrete example we may consider the parabolic 
case U — (a/2)(ri — r2)^. Another example is the Morse potential defined in the introduction. 
The pair of particles will form dumb-bell like configurations. Then the motion consists of two 
independent parts: The free motion of the center of mass R = {ri + /2 and the mass velocity 
V = (vi +V2)/2. The corresponding coordinates are Xi — (xn +a;2i)/2 and X2 — {xi2 +X22)/2. 
The relative motion under the infiuence of the forces is described by the relative radius vectors 
r = (ri — r2) and the relative velocity v = {vi — 1)2). The relative coordinates are xi = {xn — X12) 
and X2 = {xi2 — X21). Let us first study the deterministic equations. The motion of the center of 
mass is described by the equations 



iV = 



The relative motion is described by 



1 r 



F V 



F V - 



(20) 



C/'(r)- 



[F{V 



F{V -v)] 



(21) 



This system possesses two types of attractors. The first one corresponds to a translational motion 



were the two particles move nearly parallel and we have <C V'^ ~ Vq 
find in quadratic approximation in v: 

V = ["/{V^) + y{V^)v^] V - 2y{V^)iV ■ v)v 



With this assumption we 



(22) 



Assuming that the term of second order O(w^) is bounded and remains small we may conclude 
from this dynamical equation that all velocity states V"^ ~ Uq will converge to the attractor state 
= Vq. Since j'{vq) > the coupling to the relative motion may lead to an enhancement of 
the translation due to an energy fiow from the relative to the translational mode. For the relative 
velocities we get in linear approximation 



mv + U'[r)- = T ■ 
r 

Here T is a tensor which we call friction tensor 

a2 



— +-i{V^)]5 + 2i{V^){VV) 
m 



(23) 



(24) 



(25) 



From the dynamical equation for v we find for the relative energy 

1 ^2 ^ U{r)) =-vT-v + 0{v^) 

For translational velocities near to the root i.e. ~ fg the tensor T is positive, having only 
positive eigenvalues. In this case the r.h.s. of the energy equation is negative i.e. the energy tends 
to zero, both particles collapse to the minimum of the potential. In other words the attractor of 
motion is 



V = VGn] R{t) = vont + R{0) 



(26) 
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In the attractor state V'^ = Vq itself the (positive) tensor reads 

r^- + 2j\vi){VV) (27) 
m 

This means, excitations with v perpendicular to V (i.e. (VV) = 0) are only weakly damped and 
excitations with v parallel to V show a stronger damping. With decreasing < Vq the tensor 
r shows a bifurcation. This happens, when the first eigenvalue crosses zero, at this point the 
translation mode is getting unstable. As a consequence the terms O(w^) may increase unboundedly 
and the translational mode breaks down. A similar bifurcation has been found for the one- 
dimensional case in 

Let us consider now the influence of noise, naturally we expect some distribution around 
the attractors. In the stable translational regime we flnd in some approximation the stationary 
distribution 

/(°)(V,t;,r) = c(l + ^vA '° cxp 

This corresponds to a driven motion of a free particle located in the center of mass supplemented 
by a small oscillatory relative motion against the center of mass. The solutions for the rota- 
tional model are similar to what we have found for the case of external fields. The probability is 
distributed around two limit cycles corresponding to left or right rotations. 

Summarizing our finding we may state: For two interacting active particles there exists a 
translational mode in which the center of mass of the dumb-bell makes a driven Brownian motion 
similar to a free motion of the center of mass. In the rotational mode the center of the dumb-bell 
is at rest and the system is driven to rotate around the center of mass. In this mode only the 
internal degrees of freedom are excited and we observe driven rotations. The attractor region of 
the translational state has been described above, the attractor region of the rotational state and 
the exact distribution function has still to be explored. 



(28) 



5 The Dynamics of self-confined swarms of active Brownian 
particles 

The collective motion of more complex iV-particle systems (swarms) was studied already in a series 
of papers by Vicsek et al. based on a spin glass model in the velocity space 00]. First studies 
of the dynamics of harmonic swarms of driven Brownian particles were reported recently |20II27| . 
The dynamics of clusters of active particles with Morse interactions was studied first in ■ Here 
we restrict ourselves to Morse forces, generalizing our findings reported in the previous work j2f)| . 
Let us consider first the case that all particles form one cluster which is hold together by relatively 
long range Morse forces. Approximating the sum of forces in a mean field approximation we may 
represent the local field near to the center of the swarm 

by an anharmonic oscillator potential centered around R, the center of mass. In the normal 
representation the potential reads 

U{xi,X2) = ^ [ai{xi - Xi) + a2{x2 - X2)] (30) 

In arbitrary representations we have a second-rank tensor a with the eigenvalues ai,a2, which 
determine the normal modes of oscillations around the center. The tensor depends on the pa- 
rameters of the force and on the shape of the swarm. In this linear approximation the Langevin 
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equation for an individual Brownian particle in a iV-particle system bears some similarity to the 
parabolic case studied earlier |2UII27j : 

r^=Vi , mv, + a{r,- R{t)) ^ -mj{v^^)v,~ a'^{v,-V) + V2D^i{t) (31) 

The difference to the case studied before is the anharmonicity and the existence of a weak global 
coupling. Neglecting first the effects of anharmonicity we may predict, on the basis of the previous 
findings, that clusters with rotational symmetry will show a translational and a rotational mode 
very similar to the case N = 2. Therefore we expect to see translating as well as rotating clusters 
which might change their sense of rotation due to stochastic effects (see Fig. ^ . With increasing 
asymmetry of the shape, rotating clusters are expected to get unstable similar as asymmetric 
driven oscillators |2S]. In order to check these predictions at least qualitatively we have carried 
out several simulations for swarms with Morse interactions (see Figs. Handl^ll. 




Figure 1: Rotating cluster of 20 particles for different time steps. The arrows correspond to the 
velocity of the single particle. Because of the influence of noise the cluster changes the direction 
of rotation randomly from I 2(j| . 

The basic results of our observations may be summarized as follows: 

7 




cluster drift We see in tlie simulations drifting clusters rotating very slowly and clusters without 
rotations which move rather fast. This corresponds to the translational mode studied for 
N = 2. Here most of the energy is concentrated in the kinetic energy of translational 
movement. 

generation of rotations As we see from the simulations, small Morse clusters up to ~ 20 gen- 
erate left /right rotations around their center of mass. The angular momentum distribution 
is bistable. This corresponds to the rotational mode studied above for iV = 1, 2. 

breakdown of rotations The rotation of clusters may come to a stop due to several reasons. The 
first is the anharmonicity of clusters. As we have shown above |26| . strong anharmonicity 
destroys the rotational mode. Another reason are noise induced transitions, this will be 
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investigated in a forthcoming paper [28j. 

shape distribution With increasing noise the shape of the clusters is amoeba-like and is getting 
more and more complicated. A theoretical interpretations of the shape dynamics is still 
missing. 

cluster composition With increasing noise we observe a distribution of clusters of different size. 
Again a theory of clustering in the two-dimensional case is still missing. For the case of one- 
dimensional rings with Morse interactions several theoretical results are available |24II25 | . 

6 Conclusions 

We studied here the active Brownian dynamics of a finite number of confined or self-confined 
particles with velocity-dependent friction. Confinement was created 

(i) by external parabolic forces, 

(ii) by attracting Morse interactions. 

We have given here first an analysis of several simple cases as the motion of noninteracting driven 
particles in external potentials and the collective and relative motion of two driven particles. 
Based on these analytical investigations we have made several predictions for the behavior of 
Morse clusters with small particle numbers. This way we could identify several qualitative modes 
of movement. Further we have made a numerical study of special N particle Morse systems. In 
particular we investigated the rotational and translational modes of the swarm and the clustering 
phenomena. 

We did not intend here to model any particular problem of biological or social collective 
movement. We note however that the study of dynamic modes of collective movement of swarms 
may be of some importance for the understanding of many biological and social collective motions. 
To support this view we refer to the book of Okubo and Levin 29 where the modes of collective 
motions of swarms of animals are classified in a way which reminds very much the dynamical 
modes of the model investigated here. 
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